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We discuss the exact solutions of various models of the statistics of dimer coverings of a Bethe 
lattice. We reproduce the well-known exact result for noninteracting hard-core dimers by both a very 
simple geometrical argument and a general algebraic formulation for lattice statistical problems. The 
algebraic formulation enables us to discuss loop corrections for finite dimensional lattices. For the 
Bethe lattice we also obtain the exact solution when either a) the dimers interact via a short-range 
interaction or b) the underlying lattice is anisotropic. We give the exact solution for a special limit of 
dimers on a Bethe lattice in a quenched random potential in which we consider the maximal covering 
of dimers on random clusters at site occupation probability p. Surprisingly the partition function for 
"maximal coverage" on the Bethe lattice is identical to that for the statistics of branched polymers 
when the activity for a monomer unit is set equal to —p. Finally we give an exact solution for the 
number of residual vacancies when hard-core dimers are randomly deposited on a one dimensional 
lattice. 

PACS numbers: 02.10.Ox, 05.50.-|-q, 64.60.Cn 
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I. INTRODUCTION 



The statistics of covering a lattice with monomers or 
dimers has a long and continuing history in condensed 
matter physics. Recently there has been a revival of in- 
terest in this topic in connection with a number of seem- 
ingly unrelated problems such as quantum fluctuations in 
Heisenberg antiferromagnets,^ stability and dynamics in 
granular systemSfS, phase transitions in certain complex 
fluids;^ dynamics of catalysis on surfaces?^ and the bio- 
physics of membranes^ Accordingly, we have been led to 
revisit this problem with the goals of a) drastically sim- 
plifying the derivation of existing approximations and b) 
providing a framework within which the more modern 
techniques of statistical mechanics can be applied. 

The first studies of the statistics of dimer cover- 
ings of a lattice were carried out more that 50 years 
g^gQ^6j7j8 obtaining results analogous to those of the Bethe 
approximation^iiSii-l for the Ising and Heisenberg mod- 
els. At that time, the relation of this approximation 
scheme to the structure of the Cayley tree (a recursive 
"lattice," an example of which, with coordination num- 
ber q, is shown in Fig. ^| was apparently not known. 
It was later recognized by SykeSfi^ who apparently first 
coined the term "Bethe lattice," that the Bethe approx- 
imation was to be associated with local properties eval- 
uated near the center of the tree, in order to avoid sur- 
face effects which, for d-dimensional hypercubic lattices, 
are unphysical. The pathological effects of the anoma- 
lously large surface were later studied by several authors 
in the 1970'sr^ii^ii^*i& but, as was clear from the work of 
Fisher and Gaunt in 1964,1— it was the results for local 
properties at the center of the infinite tree which could 
be connected to those of hypercubic lattices in the limit 
d ~f oo. They obtained expansions in powers of 1/d for 
the coefficients of series expansions in the coupling con- 
stant. However, it was later shown by the renormaliza- 
tion groupiSi that the critical exponents for typical lattice 



models were those of mean-field theory for d > A. 

Most of the results of this paper will be obtained for 
the Bethe lattice. We implement the Bethe lattice con- 
dition either by explicitly considering sites far from the 
boundary or, alternatively, by using a formulation appro- 
priate to periodic lattices and then introducing approx- 
imations which become exact when the lattice does not 
support any loops By initially treating a periodic lat- 
tice we eliminate anomalous surface effects. We are thus 
assured that our results are characteristic of the interior 
of the tree and should be associated with what is now 
commonly called a Bethe lattice. 

Of course, an important aim is to treat real d- 
dimensional lattices. To this end there have been a 
number of papers dealing with series expansions for the 
problem of dimer or monomer-dimer coverings of a lat- 
tice. NagleSS in 1966 developed a series expansion in 
powers of the dimer activity z for a number of two- and 
three-dimensional lattices. Longer series were later ob- 
tained obtained by GaunliSi who exploited the relation 
between the dimer problem and the Ising model in a 
field. Alternative formulations, based on the generaliza- 
tion of the Mayer cluster expansionSS by Rushbrooke and 
Scoins^^ have also been giveniS4 More recently Brazhnik 
and PreedS^ have given a formulation suitable not only 
for dimers but for more complicated entities. Here we 
treat the cases of both noninteracting and interacting 
dimers. We first address these problems using a simple 
intuitive geometrical approach. Then we adopt an alge- 
braic approach based on a transformation introduced by 
Shapir^^ which enables us to develop an expansion for 
a d-dimensional lattice in which the leading term is the 
exact result for the Bethe lattice of the same coordina- 
tion number q. For noninteracting dimers corrections are 
obtained in powers of 1/q and the activity z, similar to 
the result of NagleiSS 

Briefly this paper is organized as follows. In Sec. II 
we present a simple geometrical derivation of the Bethe 
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FIG. 1: A Cayley tree with coordination number q — 3. The 
tree consists of three generations of bonds starting from the 
seed site, S. 



approximation for the statistics of distributing hard-core 
(but otherwise noninteracting) dimers on a lattice. We 
also give an alternative algebraic derivation of this re- 
sult which enables us to generate loop corrections for the 
case of a c?-dimensional periodic lattice. We also briefly 
consider noninteracting dimers on an anisotropic lattice. 
In Sec. Ill we consider a model which includes dimer- 
dimer interactions. We derive the Bethe approximation 
for this model by both the geometrical and algebraic ap- 
proaches. The exact solution of a model in the special 
limit of quenched randomness in which each cluster of 
randomly occupied sites is maximally covered by nonin- 
teracting dimers is given in Sec. IV. In Sec. V we consider 
a model of random deposition of hard-core dimers in one 
dimension. We give an exact result for the fraction of 
sites which remain vacant after deposition is completed. 
Finally, our conclusions are summarized in Sec. VI. In a 
future paper we will generalize our approach to treat the 
statistics of entities more complicated than dimers. 



II. NONINTERACTING DIMERS ON A BETHE 
LATTICE 

A. Geometrical Derivation of the Exact Result 

Here we will develop a formula for p, the average num- 
ber of dimers per edge on a lattice in terms of the dimer 
chemical potential fj,, or preferably in terms of the dimer 
"activity" z = e^^, where (3 = l/{kT). (Interactions 
between dimers and the underlying lattice as well as in- 
ternal degrees of freedom of the dimer are easily included 



by a redefinition of z.) We will give what we believe to 
be the simplest possible "geometric" derivation of the 
well-knownSiL£ "Bethe approximation," which has been 
showr.^ to be surprisingly accurate for some two dimen- 
sional lattices, and which is exact when applied to the 
Bethe lattice. We assume that only dimers can be ad- 
sorbed on the lattice, and that they can be adsorbed 
only as lying dimers (i.e., an adsorbed dimer covers two 
lattice sites). In the present section we will assume that 
there are no interactions between adsorbed dimers ex- 
cept for the hard-core restriction that two dimers cannot 
touch the same site. It is easy to generalize this treat- 
ment to deal with the two-component lattice and also to 
include the possibility of adsorption of standing dimers 
and monomers. 

As stated in the introduction, we initially deal with a 
periodic d-dimensional lattice and will obtain the results 
for the Bethe lattice by introducing an approximation 
which is exact when the lattice is tree-like (i. e. it con- 
tains no loops). As a starting point of this discussion we 
introduce the (unnormalized) probability that there are 
Nd dimers on the lattice, Wn{Nd)z^° , where Wn{Nd) 
is the number of distinct configurations of Njj dimers on 
a lattice of N sites. The equation of state for this dimer 
system is the relation between z and N o/N , where Njj is 
the average value of . In the Bethe approximatiop^^*^ 
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where q is number of sites which are nearest neighbors 
of a given site. To emphasize that this result applies to 
a Bethe lattice we express it in terms of densities such 
as p = 2N d/ (qN), the fraction of bonds (a bond is an 
edge connecting two sites) covered by a dimer or qp, the 
fraction of sites which are covered by a dimer. Stated 
alternatively, p is the probability that a given bond is 
covered by a dimer and qp is the probability that a given 
site is covered. In terms of these variables the equation 
of state is 



z = p{l-p){l-qpy 



or, equivalently. 



l + 2qz- VI + 4crz 
2{1 + q^z) 



(2) 



(3) 



where a ^ q — 1. 

For this discussion it is convenient to define some ter- 
minology to describe configurations of dimers on a lat- 
tice. In Fig. 121 we show a small section of a lattice where 
one sees occupied bonds (D) on which dimers are placed 
and unoccupied bonds, which can be of three types de- 
pending on whether or not the two sites of the bond are 
covered by dimers. If both sites of an unoccupied bond 
are uncovered, the bond is called a "dimer vacancy," (V) . 
If both sites of the unoccupied bond are covered, the 
bond is called (in anticipation of Sec. Ill where we al- 
low dimer-dimer interactions) an "interacting bond" (I), 
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FIG. 2: Classification of bonds in a lattice partially covered 
by dimers. As indicated in the text D labels a dimer, V a 
dimer vacancy, I an interacting bond, and H a half bond. 



and if only one site of the unoccupied bond is covered, 
the bond is called a "half bond" (H). We first observe 
that in the thermodynamic limit the average and most 
probable values of Nd differ negligibly. Thus we may say 
that the average number of dimers N jj is the value of No 
for which the probability Wn{Nd)z^'^ is maximal. The 
condition that this quantity be stationary is 



Wn{Nl 



I) 



1 



Wn{N 



(4) 



D 



Consider the left-hand side of this equation. From any 
configuration of No dimers we can obtain a configura- 
tion of No + 1 dimers by placing an additional dimer 
on a dimer vacancy. The number of new configurations 
obtained by placing a dimer on one of the Ny dimer 
vacancies of each configuration of No dimers is thus 
NvWn{No), where A''^ is the average number of dimer 
vacancies for configurations having No dimers. However, 
we note that in this new set of configurations each con- 
figuration occurs No + 1 times because each dimer in the 
new configuration could have been the one newly added. 
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where we will replace No + 1 by No in the thermody- 
namic limit. Applying this for the most probable value 
of No gives the simple result^ 
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If we introduce the density of dimer vacancies py by py = 
Ny/No, where No is the total number of bonds, then 
this may be written as 



Pv_ 
P 



(7) 



Note that p (or py) is defined to be the total number of 
dimer (or dimer vacancies) divided by the total number of 
bonds. These quantities can also be defined locally for a 
subsystem of bonds. It is not obvious that when p and py 
are replaced by their local versions that Eq. will still 
hold. Indeed, for the Ising model on a Cayley tree in an 
applied field, the total magnetization divided by the total 
number of sites is not the same as the local magnetization 
of a single site near the center of the tree. We now argue 
that py/p does not depend on position within the tree. 
Consider two subsystemsSS Si and S2 which consist of 
A^i and N2 bonds, respectively, where iVi and N2 are 
both large compared to 1 (so that we do not need to 
worry about discreteness effects). Suppose py/p assumes 
different values within these two subsystems. Then, if we 
move a dimer from subsystem Si to subsystem S2, the 
ratio of the number of configurations W before moving 
the dimer to that, W , after moving the dimer is 
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where the subscripts identify the subsystem in question. 
We may view the pyn for subsystem S„ as being order 
parameters. In order that Wn{No) actually be maximal 
with respect to variation of these order parameters, it 
must be stationary with respect to moving a dimer. So 
W/W = 1, which implies that py/p is the same for all 
large subsystems. Thus Eq. Q, which initially involved 
the ratio of global properties, is actually valid when in- 
terpreted as a relation between local quantities. 

We now express py as an explicit function of p in a 
region far from the boundary, so that Eq. Q yields 
the equation of state we seek. If A and B are nearest 
neighboring sites, then 

py = P (sites A and B are vacant) 
= P(A is vacant) 

xP(B is vacant, given that A is vacant) 

= (i-qp) 

xP(B is vacant, given that A is vacant) . (9) 

Here the symbol P() denotes the probability of the event 
inside the parentheses. The basic probability space con- 
sists of all configurations of dimers on the lattice. The 
(unnormalized) probability assigned to a configuration 
is , where A'^^i is the number of dimers in the con- 
figuration. Equivalently, we could limit our probability 
space to the set of configurations containing exactly No 
dimers, with all configurations equiprobable. The condi- 
tional probability [the last factor in the last line of Eq. 
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is denoted P(B vacant | A vacant) and is, by defini- 
tion, equal to P(B vacant and A vacant)/P(A vacant). 
More generally, for any two events Ei and E2 



P{E2\Ei) = P{E2 and Ei)/P{Ei) 



(10) 



We now invoke the approximation in which we re- 
place the factor P(B vacant | A vacant) in Eq. ((HJ by 
P(B vacant | bond AB unoccupied). This replacement 
ignores the possibility that at most one of the other bonds 
touching site A might be occupied and thereby might 
indirectly affect whether site B is occupied or not. Of 
course this possible inaccuracy can only arise if there is 
some indirect path from site A to site B not going through 
the bond AB. Since this approximation (which we will re- 
fer to as the tree decoupling) is exact for the Bethe lattice 
it leads to the Bethe approximation. Thus 

py = (1 - qp) 

X P(B is vacant | bond AB is unoccupied) (11) 

The probability that a particular bond is unoccupied is 
1 — p, so that 

P(B vacant | bond AB unoccupied) 
P(B vacant and bond AB unoccupied) 

P(bond AB unoccupied) 
P(B vacant and bond AB unoccupied) 



(12) 



If B is vacant, the bond AB must be unoccupied; there- 
fore the numerator of Eq. is just P(B vacant) = 
1 — qp. Thus 



Pv 



(l-gp)V(l-p) 



(13) 



which, in combination with Eq. (Q, leads to 
Eq. We have not succeeded in construct- 

ing a simple argument to estimate the magnitude of 
the difference between P(B vacant | A vacant) and 
P(B vacant | bond AB unoccupied) for d-dimensional 
lattices. Qualitatively, the important point is that, 
within the family of configurations in which bond AB 
is unoccupied, the probability that site B be vacant is 
not significantly influenced by the presence of a dimer on 
(at most) one of the other bonds emanating from A. 

In Appendix C we present a simple calculation of the 
major correction to the tree approximation on a planar 
triangular lattice. However, the calculation is not the 
first term in a systematic series and is not easily extended 
to the square lattice. In the next section we will present 
a formalism which enables one to systematically generate 
corrections to the tree approximation. 



B. Solution by Construction of an Effective 
Hamiltonian 



We now apply a technique introduced previously-'J- in 
order to treat here the nonthermal statistical problem of 



constructing a generating function for covering a lattice 
with 0,1,2,... hard-core dimers. We wish to identify this 
generating function with a partition function, Z of the 
form Z — Tt exp(— /37i), where Ti. can then be interpreted 
as the Hamiltonian for the statistical problem. 

The first step in this program is obviously to construct 
the effective Hamiltonian. This can be done by writing2& 



(14) 



where (ij) indicates that the sum is over pairs of nearest 
neighbors. For the expansion of the partition function 
in powers of z to count all possible dimer configurations, 
where z is the dimer activity, the following trace rules 
are imposed on the operators sf. 



Trj;s" — Cn 



n = 0,1,, 



(15) 



where Tr.^ indicates a trace over states of site i and we 
set Co = Ci = 1 and C„ = for n > 1. It is not ac- 
tually necessary to explicitly construct such an operator 
because the only property of these operators we need in 
order to construct the partition function is the trace rules 
of Eq. H15|l . We see that the fact that the trace of two 
or more operators at the same site vanishes, implements 
exactly the hard-core constraint for dimers. Thus the 
partition function Z = Ti exp(— /37i) will indeed give the 
grand partition function for dimers as a function of their 
chemical potential fi: 



(16) 



where the sum is over all configurations C of dimers and 
n{C) is the number of dimers present in the configuration 
C. From Z we can get the fraction of sites covered by 
dimers p as a function of the dimer chemical potential 
via 



2Nd _ 2 d\nZ _ 2 d\nZ 
Nq ~ Nqj3 dp, ~ Nq dlnz 



(17) 



Since the "spin operators" Si commute with one another, 
we have a mapping of the athermal problem of dimers on 
a lattice into a statistical mechanical problem involving 
classical spins with a given Hamiltonian. 

Now we use this mapping to a) construct the exact so- 
lution for the partition function for a Bethe lattice and b) 
generate series expansions for finite dimensional lattices. 
To do that we develop a perturbation theory for a peri- 
odic lattice in which the leading term contains the sum 
of all contributions from tree diagrams. For this purpose 
we write Z = TiY[(ij) fijj where = l + zSiSj [the trace 
rules allow us to linearize the exponential in Eq. (|14|) ] or 
equivalently 



Tr<! 



n— 



(18) 
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where g can be chosen arbitrarily. We evaluate this per- 
turbatively as 



Thus 



Z = Tr- 



(ij) 



(19) 



where Vij — fij/igigj) ~ 1- We expand in powers of A 
which we set equal to unity at the end. Each term in 
this expansion which involves at least one power of A can 
be associated with a diagram in which the factor Vij is 
associated with a line connecting sites i and j. We now 
choose g so that diagrams having at least one line which is 
connected to only a single site (j) give zero contribution. 
Since all diagrams on a Bethe lattice have at least one free 
end, this choice of g will lead to an exact evaluation of 
the partition function for a Bethe lattice and will enable 
us to generate loop corrections for d-dimensional lattices. 
The condition we implement is that 



Tr, {g]V,,) = 
which can be written in the form 



9i 



Tr 



(20) 



(21) 



where a = <? — 1. This is a nonlinear equation for the 
function gj, but in view of the trace rules it is easily 
solved. From the form of one sees that gi has to be 
of the form 



= A + Bsi . 
By substituting this form into Eq. I|21|l : 

Tr,- [{l + zs,Sj){A + BsjY] 



A + Bsi = 



Tvj{A + Bsj 



Using the trace rules we rewrite the right-hand 
this equation so that 



A + Bsi = 



A"" +aA''-^B + sAzA'^) 



Ai + qA'^B 
This gives rise to the two equations 

A'' + aA''-^B A + aB 



A = 
B = 



Ai + qA^B A^ + qAB 
zA" z 



At + qA-^B A + qB 
We may solve Eq. H25f) for B as 

B = ^'^^ 
a - qA^ 

Substituting this into Eq. (|26|l leads to 
A^ (l + (^z) - A^ {2qcrz + 1) + zcr 



(22) 

(23) 
side of 

(24) 

(25) 
(26) 

(27) 

(28) 



A' = 



2qza + 1 + \/r+4zCT 
2(l + g2z) 



(29) 



(We chose the positive sign before the square root to 
ensure that A ^ 1 asz^O.) Then for the Bethe lattice 
we have the exact result 



Zi/^ = Tr,g^ = TT,[A 
= A« + qBA'^ . 



Bs 



(30) 



After some algebra it can be shown that when this result 
is inserted into Eq. (|17|l we recover the result of Eq. ^ . 

Unfortunately, this formalism leads to an expansion 
of the partition function Z , whereas for d dimensional 
lattices, we would prefer to have an expansion for the 
free energy per site, F = {l/N)\n.Z . For that purpose 
we consider an expansion of the quantity Z", which in 
the limit n^Oisl + nlnZ=l + nNF. To obtain 
we introduce the n-replicated Hamiltonian 



nnn 

a=l (ij) 



(31) 



where Sia and Sip are independent operators for a 7^ /3 
and for each replica index Sia obeys the same trace rules 
as in Eq. (|15|) . Usually replicas are introduced to perform 
the quenched average (over In Z) for random problems in 
which case the averaging leads to interactions between 
different replicas. Here we introduce replicas simply to 
facilitate construction of an expansion of In Z and there 
are no interactions between different replicas. The parti- 
tion function, Z-^c-p associated with the replicated Hamil- 
tonian is 



(32) 



Thus we solve Eq. with = na(l + ^SiaSja)- 

Because different replicas are independent of one another, 
the solution to Eq. H21() is of the form 



Ds,, 



Cg^ 



(33) 



Because we need the partition function Zicp to order n^^i 
we must evaluate C up to linear order in n but D can be 
evaluated for n = because it always appears in connec- 
tion with a sum over replica indices which give a factor 
of n. The terms in gi in Eq. (|21|l independent of give 



C = TrjC"".gJ/Tr,C'«g 



or 



- Tr,,gJ/Tr,-5: 
= [l + aDriil + qDY 
which to linear order in n gives 



^ n/l + aD 

C = 1 + - In 

2 \l + qD 



(34) 



(35) 



(36) 
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The terms in gi in Eq. H21() linear in Sia give 



CD 



For n ~ this is 



Si) 



(37) 








D = z{l + aD)-^ 



(38) 



which gives 



D = 



-1 + r 
2cr 



(39) 



where r = VI + 4crz. 

Then the first term in the expansion of F is 



F 



d 

dn 



l + qD 



n=0 

+ hl(l + qD) 



(40) 



After some algebra one can show that when this is sub- 
stituted into Eq. (|17|l we recover Eq. (2). 



C. Loop Corrections 

Here we consider the expansion in powers of Vij = 
fijl{9i9j) ~ 1- Our first objective is to show that contri- 
butions to Z^e-p from disconnected diagrams are of order 
or higher and hence can be dropped. To see this we 
write 



c- 




-1 



(41) 



The contribution to Z-^cp from a disconnected diagram is 
simply the product of the contributions from each con- 
nected component. We now argue that the contribution 
to Z-cep from a single connected diagram is of order n. 
Note that Vij is a multinomial in the Sq's whose con- 
stant term is proportional to n. Thus the contribution 
to Zj-cp will get at least one factor of n, either from the 
constant term in a Vij , or from a sum over replica in- 
dices from terms in a Vij involving an s operator. Thus 
using the rt — > limit of the replica formalism we have 
eliminated unlinked diagrams. 

Now we consider the leading loop corrections to Z". 
These come from the smallest loops that can be drawn 
on the lattice. In Fig. O we show these loops for both 
the d dimensional hypercubic lattice and the plane tri- 
angular lattice. Here we will explicitly evaluate only the 
corrections from the smallest possible loop diagrams. 

We will now show that this expansion involves evalu- 
ating a modified dimer partition function for diagrams 
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FIG. 3: Leading loop corrections to for the hypercubic 
lattice (top) and for the plane triangular lattice (bottom). 



with no free ends. The contribution to the replica parti- 
tion function from such a diagram F is given by 



Kiev {.y)er ^^^^ 
= n5F{T) + 0{n^) . 



(42) 



The notation i e F means that the site i is a site covered 
by at least one bond of F. We do not change the result 
for 5F{T) if we divide this by a quantity which differs 
from unity by terms of order n. So, for later convenience 
we write 



Tr a.r^m 



l(y)Gr 



(5Zrop(F) 



Trn.erSf 



.(43) 



We now expand the product over bonds, into its 2^^'^'"^ 
terms, where 7Vb(F) is the number of bonds in the set F. 
In so doing note that when considered as a multinomial 
series in {sia}, each term has the limiting value unity as 
n ^ (but the sum of all 2^^ terms is zero in this limit) . 
Since we are interested in the limit n 0, we consider 
d6Z,-cp{T)/dn (evaluated at n = 0) and write 



rf^^rep(r) 

dn 



■yer 



(_l)A'B(r)-w«(7)G'(^) ^ (44) 



where the sum is over the 2^^*^^") — 1 nonempty subsets 
7 of F, including 7 = F. (The term corresponding to the 
empty set is unity and therefore drops out when differ- 
entiated with respect to n.) Here 



G(7) = 



dn 



Ti-n 



iG7 ■ 



(45) 
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Because we divided by the factor in the denominator, 
it is no longer necessary to involve sites in F which are 
not in 7. The subtractions of subdiagrams indicated in 
Eq. (|44|) defines the cumulant operation (indicated by a 
subscript "c"), so we write 



6F{T) = Gc(r) 



(46) 



where Gc{T) is the right-hand side of Eq. This 
representation is not very efficient because it contains 
2^b{t) _ ^ terms, many of which are either zero or are 
identical to one another. The following equivalent recur- 
sive definition is more convenient: 



G,(7) = G(7)-5]'g,(7'), 



(47) 
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where the prime indicates that in the sum over subsets 7' 
we do not include 7' = 7. For the smallest loop on a hy- 
percubic lattice, namely a square of four bonds, there are 
no nonzero subtractions and Gc(7) is equal to its "bare" 
value G(7). More generally Eq. H47(l has many fewer 
terms than Eq. (|44(l and furthermore, all the cumulants 
of the subdiagrams will have been previously calculated 
in a lower order calculation. So, to implement the cumu- 
lant subtraction we only need to subtract the cumulant 
contributions of subgraphs with no free ends. 

It remains to discuss the calculation of G(7). We take 



and 

Then 
Gil) 



dn 

d_ 

dn 



C-2Ar,(7)n„Trn.e.5Ln 



(48) 



(49) 
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n«Trn. 



^-2NBh)(^l + gi:>)-^=(^)"Q(7)" 



(50) 



where Ns{j) is the number of sites in 7 and Q{"f) is a 
partition function for the graph 7: 



Qh) 



Tr 



Tr 



Y[ii+Ds.r n 



1 + ZSiSj 



{l + Dsi){l + Dsj 



ZS^Sj j 



,(51) 



where qi (7) is the number of sites neighboring to i which 
are connected to i by a bond in 7. Because of the trace 
rules, the product over sites mimics a site-dependent 
monomer activity Zi = 1 + [q — qi{j)]D and we there- 
fore have 



Qii) 



n 



Q{l\z.^ = z/{z.,Zj)) 



(52) 



where Q is the grand partition function for the set of 
bonds 7 in which the bond {ij) has the bond-dependent 
activity zij = z/{ziZj): 



Q(7;{%}) 



Tr Jl {l + z,jSiSj) 



(53) 



Thus the rcnormalized free energy associated with a dia- 
gram with no free ends is only slightly more complicated 
than its unrenormalized [gi = 1) value. To summarize: 



G(7) = -Ns{l)Hl + qD) + NB{l)\n 



1 + qD 



J2^nil + \q-q,ij)]D) 

iS7 



1 + 

lnQ(7;{z,,}) (54) 



It is a remarkable fact that for a diagram 7 with no loops, 
G(7) vanishes and this forms a nice check of computer 
programs used to evaluate G(7) for an arbitrary diagram. 
(This is easy to check for small diagrams.) Furthermore, 
for a diagram with loops and which has a free end, G(7) 
does not vanish, but its cumulant Gc(7) does vanish. 
(This is also a nice check of computer programs.) 

For the hypercubic lattice we consider the leading cor- 
rection from a square of four nearest neighbor bonds, 7. 
So we use Eq. (|54|l with q = A and (7^(7) = 2, so that 



+ 2Df. Then ^(7) = 1 + 4^ 



2z?- and 

'-J 



G(7) = 41n(l + 2D) -41n(l + 4D) + 41n 



1 + 4D 
1 + 3L» 



-In 



1 



Az 



2z' 



{\ + 2DY (l + 2£i)4 

= -41n(l + 3D) 

+ In [(1 + 2Df + 4z(l + 2Df + 2z'^ 

Now use l + 'W = z/D from Eq. so that 

Gil) 



(55) 



In 



1 



In 



z 

8 



4— (1 

z 



z 



D 

1 + — 

z^ 



(56) 



For z small D is proportional to z and this diagram gives 
a contribution to the free energy of order z*. Since there 
are d(d — l)/2 squares per site, the perturbative contri- 
bution to the free energy per site is 



5F = id(d - 1) In 



1 + - 



(57) 



The dimer density then follows using Eq. H17|l and the 
results are given in Table I for a square lattice. 

We now identify the expansion parameters in this for- 
mulation. It is clear that the free energy, F, is obtained 
as a sum of contributions associated with diagrams hav- 
ing no free ends, the smallest of which are shown in Fig. 
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TABLE I: First Loop Correction for the Square Lattice for 
the density of dimers (tt = No/N) as a function of z. Here 
TTo is the value for the Bethe lattice for q = 4, tti is the value 
when the first loop correction, Eq. 15711 . is included, and tte 
is the "exact" result obtained by extrapolation- of the series 
of Gaunt 1^ 



z 




To 




TTl 






U.UUOloU 


0.010000 


0.010000 


0.010000 


n m 07/19 

U.UIU / 





UzUUUU 





UzUUUU 





UzUUUU 


U.UZOlOD 





039999 





039999 





r\ A r\r\r\r\ 

040000 


U.UO (0/0 





059997 





UdUUUU 





r\cir\r\r\r\ 

UdUUUU 


U.UO'4:'4rlO 





0(999U 





UoUUUU 





UoUUUU 


n 07/1 1 QPi 

u.u ( ^lyo 





nnnn'70 
0999/8 





lUUUUU 





lUUUUU 


u.uy / oo 1 





119954 





119999 





120000 


U. 1Z040Z 





139917 





139998 





14UUUU 


u. looy 10 





loyoDz 





loyyyD 





IdUUUU 







179784 





179991 





180000 


n 9/1Q1 1 1 





199679 





199982 





200000 


u.oiuoy i 





219543 





219968 





220UOU 


U.oo / / DD 





2o93/U 





2^9945 





24UUUU 


U-'^iODl 10 





259158 





259909 





260000 


U.DlO / Zo 





278904 





279857 





280000 


0.782881 





298607 





299783 





300000 


1.012941 





318271 





319684 





320000 


1.335914 





337903 





339556 





340000 


1.807776 





357516 





359397 





360000 


2.533661 





377128 





379204 





380000 


3.730072 





396761 





398970 





400000 


5.902289 





416436 





418682 





420000 


10.466788 





436180 





438318 





440000 


22.802101 





456068 





457892 





460000 


81.673705 





476432 





477641 





480000 



13 This development will lead to an evaluation of F as a 
power series in the activity: 



F 



(58) 



From Eq. 14()|l we see that for the Bethe lattice and for 
large q Fn ^ /n. We assert that the contributions 
to Fn from a diagram 7 [which we denote (5F(7)] are of 
order 



5Fn{l)lFn ^ q 



-r(7) 



(59) 



where r(7) is an integer which increases with the size and 
complexity of the diagram. For instance, for a square 
of four bonds Eq. (|57|l indicates that r = 2 and for a 
loop of 2n bonds a similar result shows that r ~ n. So 
our diagrammatic formulation generates corrections in 
inverse powers of q. Furthermore, if one expands and 
gi in powers of z, one sees that Vij is of order z. This 
means that a diagram 7 with Nb{'~i) bonds contributes 



TABLE IL First Loop Correction for the Triangular Lattice. 
The notation is as in Table I, but where vri includes the loop 
correction of Eq. I6OII and ttg includes the loop correction 
derived geometrically in Appendix C. 
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0.009999 


0.009999 


0.010000 


0.003459 





011206 
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029999 
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029999 
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043360 





090121 





089991 





090000 





089987 





066712 





120261 
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119962 





097473 





150457 





149943 





150000 





149912 





138811 





180705 





179892 





180000 





179831 





195705 





210994 





209822 





210000 





209712 





276329 





241312 





239741 





240000 





239559 





394826 





271650 





269667 





270000 





269386 





577329 





302001 





299629 





300000 





299221 





876264 





332353 





329663 





330000 





329101 


1 


408720 





362668 





359786 





360000 





359056 


2 


477977 





392849 





389972 





390000 





389084 


5 


067826 





422724 





420118 





420000 





419135 


13.932546 





452078 





450085 





450000 





449161 


91.850988 





480860 





479893 





480000 





479338 



to F{ci) at order and higher. So our formulation 

involves the two expansion parameters z and \/q. 

An entirely analogous calculation gives the leading 
loop correction for the triangular lattice from triangles 
as 



SF = Nt In 



1 - 



r3 



(60) 



where Nt = 2 is the number of triangles per site. Thus 
we see that the leading correction to the Bethe lattice re- 
sult has the opposite sign for triangular lattices as com- 
pared to hypercubic lattices. The results based on Eq. 
are given in Table II. 
The present development appears to be related to that 
of Ref. '23, but the detailed relationship of the two ap- 
proaches is not clear to us. We note that the disconnected 
diagram of Fig. 3d of Nagle2& (which "is required" for 
his dimer series) does not appear in our approach. 



D. ANISOTROPY 

We may generalize the above model to allow for dif- 
ferent activities along different coordinate axes. For that 
purpose we relate the square lattice to a Bethe lattice 
in which each site is surrounded by four bonds, two of 
which we arbitrarily label as "horizontal" (or x) bonds 
and the others as "vertical" (or y) bonds. Dimers along 
the a axis have activity Za for a = x,y. The fraction of 
a bonds covered by dimers will be denoted pa and the 
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fraction of a bonds which are dimer vacancies will be 
denoted pva 



1. Geometrical Approach 
The geometrical reasoning used before yields 

•2a ^ = PVa/pOi ■ (61) 

The tree approximation, which is exact for the Bethe 
lattice, then yields 



PVcc = W/{1 - Pa) , 



(62) 



where W^^"^ is the probability that a given site is vacant, 
where 



(63) 



so that 



Thus 



and 



Px(l - Px) _ Pyi'^ - Pv) 



= W . (64) 



Pa = [1 - Vl - 4^«H^]/2 



(65) 



W = Vl - ^^xW + VI - ^ZyW - 1 . (66) 



This equation allows one to calculate as a function of 
Zx and Zy, from which all the other relevant quantities 
can be obtained. 



2. Algebraic Approach 
We use the Hamiltonian 



-I3H 



n [1 



n [1 



(67) 



where the sum "G H" means wc sum over horizontal (x) 
bonds and similarly for "e V." The operators Sj obey 
the same trace rules as before. We write the partition 
function as 

Z = T^^Ualdiylli^ + Vij], (68) 

i {ij) 

where for horizontal bonds 

Vij = [1 + ZxSiSj]/gixgjx - 1 (69) 
and for vertical bonds 

Vij = [1 + ZySiSjj/giygjy - 1 . (70) 



Now wc expand in powers of Vij and require that di- 
agrams with either a vertical or a horizontal free end 
vanish. For horizontal bonds we require that 



T^^jdjxgjyVij = , 



so that 



_ T^T^jgjxg]y[i + Z,rS^Sj] 
9ix — rr 2 2 

For vertical bonds we require that 



so that 



'^39%9%Vij = , 



'^j9'jx9jy[^ + hjSjSj] 

~ Tr o2 o2 

^^j9jx9jy 

These equations have a solution of the form 

9ia ~ ~t~ ^ot^i • 

Thus 



(71) 



(72) 



(73) 



(74) 



(75) 



-^x ^x 



Tlj[l + Za;SiSj][Ax + BxSj][Ay + BySj 

TrjlAx + BxSj]'^[Ay + BySj]'^ 



AxAy + BxAy + 2AxAyBy + SiZxAxAy 

AlAl + 2AxBxAl + 2AlAyBy ^ ^ 



and 



Trj[l + ZySiSj][Ax + BxSjf[Ay + BySj] 



i:Tj[Ax + BxSjY[Ay + B. 



Jyi=-j\ 

_ ^X^V + '^^xBxAy + A^By + SiZyA^Ay 

■^x-Ay "I" "^AxBxAy + 2A^AyBy 

Thus if we set Bx = bxAx and By = Ayby, then 

^2 ^ l + bx + 2by 
l + 2bx + 2by 

^2 _ 1 + 26^ + by 



(77) 



Aibx = 



Thus 



6„ = 



1 + 26^ + 2by 



1 + 26^ + 2by 



l + 2bx + 2by 



l + bx + 2by 
fj| 

l+2bx + by 



(78) 



(79) 



These equations reproduce those from the geometrical 
approach if one makes the identification 



Pa — Oa 

/[l + 2bx+2by] 



(80) 
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and 

W ={l + 2b^ + 2&y)-2 . 
The density of dimer bonds along a is given by 
1 dZ 

= z^-^\nlAlAl[l + 2b^ + 2by]] 



(81) 



(82) 



This can be shown to be equivalent to Eq. (|64|l . Thus 
we conclude that the algebraic approach agrees with the 
much simpler geometrical approach. But, in principle, 
the algebraic approach can be used to generate correc- 
tions to the tree approximation for d dimensional lattices. 



III. INTERACTING DIMERS 
A. Geometrical Derivation of The Exact Result 

Here we consider the case when two dimers separated 
by a single bond (an interacting bond as shown in Fig. 
1) have an interaction energy —a. Thus the energy of 
a configuration of dimers on the lattice is —Nja, where 
Nj is the number of interacting bonds. We start by ex- 
pressing Ni in terms of Ny- For that purpose we record 
the following sum rules for periodic lattices. The first 
sum rule expresses the fact that each bond is uniquely a 
member of one of the four sets shown in Fig. 1, so that 
the total number of bonds of the lattice Nb is given by 



Ni + Nd + N 



Nf 



Nq/2 



(83) 



where Nh is the number of half bonds (see Fig. 1). The 
second sum rule is obtained by imagining putting a stick 
on each of the q bonds emanating from each occupied 
site. The total number of sticks is obviously q{2ND)- 
Each half bond has one stick, whereas each interacting 
bond and each dimer bond is covered by two sticks, so 
that 

2qND ^ Nh + 2Ni + 2Nd ■ (84) 

Combining these two relations we get 

Ni ^ -~Nq/2 + {2q - 1)Nd + Nv ■ (85) 

Thus the probability P{Nd) that there are Nu dimers 
on the lattice is 

P{Nd) - Q-i^e'='"^^z^° 
Q 

= const X [ze^"(2,-l)]Ar„^g/3aiV.(e) ^(gg) 

9 

where Q is the partition function, the sum over Q runs 
over all configurations with Nd dimers on the lattice, and 



Nv{Q) is the number of dimer vacancies in the graph Q. 
Similarly 

P{Nd + 1) = const. X [ze''"(29-l)]JVz, + l ^g/3aArv(e') ^ 

Q' 

(87) 

where the sum over Q' runs over all the graphs with 
Ne) -t- 1 dimers on the lattice. Let G"{G) be a graph 
formed by adding an additional dimer to a graph Q. The 
additional dimer is, of course, placed on one of the dimer 
vacancies AB in the graph Q. The number of vertices 
adjacent to A (excluding B) is cr, and the number of ver- 
tices adjacent to B (excluding A) is also a. If m of these 
2(7 vertices are vacant, then the number of dimer vacan- 
cies in the graph Q" {Q) is Nv{Q) — m — 1. If we look at 
all the different graphs G"{G) which can be generated by 
adding a dimer to a particular graph Q, m assumes the 
values (1, .., 2(t) with the respective probabilities p{m) 
which we shall calculate. Furthermore, since it is quite 
clear that the configurations of two small subsections of 
a large graph Q which are distant from each other are 
statistically independent, we assert that the probability 
distribution p(m) is the same for almost all graphs S (but 
p{m) does depend on Nd/N which has the same value 
for all the graphs Q). Thus for almost all graphs G we 
have 

^ ^PaNv[g"{g)] 
e"(e) 



2<7 



A^v(e)e'3"^-(^) ^p(m)e 



(88) 



If {G"{G)} is the set of all graphs which can be gener- 
ated by adding a dimer to a particular graph G, and 
J2g{S"{G)} is the set which is the union of all the 
sets {G"{G)}, then J2g{S" {G)} is identical to the set of 
graphs {G'}, except that each graph G' occurs Njj + 1 
times in the union (since each graph G' has N^ + 1 " an- 
cestors" G which can be obtained by removing a single 
dimer). Thus we obtain 



{Nd + 1) 5Z 6' 



g g"(g) 

( |^p(m)e-'3"(™+iM ^7Vy(^)e^"^-(^) 

\m=0 / g 

p(m)e-'3"('"+i)^]Vv.(Ari3) e'^"^-^^) 



Vm=0 



(89) 



where 



NviNo) = ^A^y(e)e^"^^(^V^e^"^^(^) .(90) 

a g 
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FIG. 4: Left: the bond AB and its environment in a Bethe 
lattice with q — 4. The subsystems Ta and Tb are the regions 
outside the dotted hues. We show the sites neighboring to A 
and B, which are numbered 1 to 6 (1 to 2a). There are no in- 
direct paths connecting any pairs of numbered sites which do 
not pass through A and/or B. Middle: Environment of a bond 
AB on a square lattice where we show indirect paths which 
connect adjacent numbered sites, in which case the strong 
form of tree decoupling is not exact. Right: Environment 
of a bond AB on a triangular lattice. Here the decoupling 
is worse than for a square lattice because a) there are near- 
est neighbor bonds connecting numbered sites and b) some 
neighbors of A are simultaneously neighbors of B. 



For a given value of z, the value of Njy which maximizes 
P{Nd) satisfies the condition P{Nd + 1)/P{Nd) = 1, 



1 = -.M2,-1) N_v{^d) 



la 

E 

m=0 



p{m)e 



-i3a{m+l) 



(91) 



(Note that when a = this reduces to Eq. © for nonin- 
teracting dimers.) The preceding statements are true on 
all lattices with q nearest neighbors. 

We now obtain Ny and p{m) explicitly as functions 
of Njj, so that Eq. H91|l becomes an explicit equation 
of state for dimers. We do this using an approximation, 
which for reasons explained below we call the "strong 
form of tree decoupling," in which we replace the d- 
dimensional lattice by a Bethe lattice. To obtain p{m) 
we first relate p{m) to the quantity Pcon which we de- 
fine to be the conditional probability that site i is vacant 
when site j is known to be vacant, where sites i and 
j are nearest neighboring sites. Recall that A and B 
are vacant neighboring sites. We label the vertices ad- 
jacent to A (excluding B) by the index i [i ~ l,..,cr) 
and the vertices adjacent to B (excluding A) by the val- 
ues i = cr -I- 1, . . . 2cr, as shown in Fig. 0] We define a 
random variable Xi which has the value 1 if the vertex 
i is vacant and the value if a dimer is touching ver- 
tex i. On the Bethe lattice the (i = 1, 2, . . . 2a) are 
independent random variables. Note that P{Xi = 1) is 
the conditional probability that a neighbor of B is vacant 
given that not only B but also A is vacant. On a Bethe 
lattice the condition that A is also vacant is irrelevant in 
this context, in which case P{Xi = 1) is just Pcon, so that 



p{m) 



P{Xi + X2 . 
m!(2fT — m) 



X2a = m) 



■Pcon(l Peon) 



2a — m 



(92) 



2a 



J2 p{m)e~'"'^^^+^^ 

(3a 



m=0 

e 



- "-""[Pcone-'^^ + l-pcon]' 



Thus, Eq. l|9T|l becomes 
Nv 



l = z=^[e + l-pconeV^ , 

1\ D 



(93) 



(94) 



where e — e^" — 1 and we replaced A^/j + 1 by iVc in 
the thermodynamic limit. Now we use Eq. which we 
write in the form 



Pv = (1 - qp)Pcoi 
Then Eq. becomes 

1 = ZPcon [e+l-pc 



l2cr 



(95) 



(96) 



It remains to determine Pcon in terms of p. 

In order to calculate Pcon (as a function of p), we first 
calculate a simpler quantity p' , the conditional probabil- 
ity that B is vacant, given that bond AB is unoccupied 
(i.e. not covered by a dimer). Note that an unoccupied 
bond is not necessarily a dimer vacancy, since a dimer 
may be touching one or both ends of the bond. We write 

p' — P(B vacant | bond AB unoccupied) 

P(B vacant and bond AB unoccupied) 
P(bond AB unoccupied) 
= P(B vacant)/P(bond AB unoccupied) 
= il-qp)/{l-p) . (97) 

On the other hand, we will calculate p' in terms of the 
weights (unnormalized probabilities) of certain graphs, 
yielding the information which is needed to calculate 

Peon • 

We focus our attention on a bond AB and we continue 
to consider the Bethe lattice. For this discussion we di- 
vide the lattice into two parts by cutting the bond AB, 
so that one part, Ta, contains site A and all sites ac- 
cessible to site A without going through site B and the 
other part, Tb, is defined similarly, as shown in Fig. 01 
On the Bethe lattice, when the bond AB is unoccupied, 
these two parts Ta and Tb are independent subsystems. 
In that case we see, from Eq. that the "weight" 

of a configuration of either Ta or Tb is obtained as the 
product of (a) a factor ( = z exp[(3a{2q — 1)] for every 
bond covered by a dimer— (the actual value of C will not 
appear in the final results) and (b) a factor e^" for every 
bond which is a dimer vacancy. We define 

wa,o — sum of weights of all graphs on T^ which 

contain no dimers touching A 
wa,i = sum of weights of all graphs on Ta which 
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Wb,0 



wba 



contain one dimer touching A 
sum of weights of all graphs on Tg which 

contain no dimers touching B 
sum of weights of all graphs on Tg which 

contain one dimer touching B . (98) 



Note that we are only concerned with graphs in which 
the bond AB is fixed (to be unoccupied) , so that weights 
of graphs on the subtrees T\ and Tb are well defined. 
Then 



P(B vacant) 
P(AB unoccupied) 



(99) 



WA,oWB,oe^°' + WA,0WB,1 + Wa,iWb,1 + WA,lWB,0 

The first term in the numerator is the total weight of 
graphs in which A and B are vacant (so that bond 
AB is a dimer vacancy), and the second term is the 
total weight of graphs in which B is vacant and A is 
occupied (so that bond AB is not a dimer vacancy). 
The terms in the denominator are similarly interpreted. 
Since P(B vacant) = P(A vacant) far from the bound- 
ary, we have wb,owa,i = wb,iwa.o- Thus wa,i/'wa,o — 
wb,i/wb,q- 

Dividing the numerator and denominator of Eq. I|99l) 
by wa,oWb,o we obtain the following relation between p' 
and u = WA^i/wAfi- 



SO that 



2u ' 



f - 
I- P 



g/3a + u2 _|_ 2u 

Now we relate Pcon to u: 

P(A and B vacant) 



(100) 



(101) 



P(A vacant) 

WAfiWBfie 



Pa 



WAfiWBfiel^" + WA,0WB,1 

e+l 



1 



(102) 



We now solve Eq. H102I) for u and substitute the result 
into Eq. HlOlfl to get 



(2e + A + l) + (e + l)/peo„ = 



where A — pcr/{l — qp). (Note that as a ^ 
becomes equal to p' , as expected.) Thus we get 



1 



qp 



PconCr 



P [ePcon - e - 1] bcon - 1] 

so that Eq. (|9(j|l can be written as 



1 



1 



-[e+l -Pconf] 



(103) 

0, Peon 
(104) 

(105) 
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FIG. 5: Schematic graphs of G(p) for e < Ec (left) and for e > 
Ec (right). In the right panel we show the line at — ln[(T 20(e)] 
at which the net area between the line and the curve of G{p) 
is zero. This is the "equal area" construction which gives 
zo{e). We also show lines for a value of z {z<) which is less 
than zo(e) and a value of z (z>) which is larger than zo[e). 
The values of pj are the values of p at which the line for zo{e) 
intersects G{p). 



ln((Tz) = G(p) 



(106) 



where 



G{p) = 21npcon(p) - hl[l - J3con(p)] 

+ (2(7 - 1) ln[e + 1 - epcon(p)] , (107) 

where pcon (p) is determined by Eq. H103|l 

The physical range of p is [0, 1/q]- As p ^ 0, Pcon 1 
and G{p) +oo. As p ^ 1/q, Pcon and G{p) — oo. 
One can show, either analytically or by making computer 
plots, that if e is less than a certain critical value Ec (which 
we shall calculate) then G{p) is a monotone decreasing 
function of p in the interval [0, 1/g]. In this case there 
is a unique solution of Eq. (|106|l for p for each positive 
value of the activity z. If e > ec, then dG/dp is negative 
for small p and for p near 1/q, but there is an interme- 
diate region in which dG/dp > 0, as is illustrated in Fig. 
O In this case there is a certain range of z in which 
there are three values of p which satisfy Eq. (|106|l . As 
explained in Appendix^ for z less than a certain value 
zo(e), which is determined by an equal- area construction, 
the physical solution (i.e. the most probable state) is a 
homogeneous state whose density of dimers is given by 
the smallest solution of Eq. (|106|l . For z > zo{e) the 
physical solution is a homogeneous state whose density 
of dimers is given by the largest solution of Eq. (|106|1 . 
For z — zo{e) the physical solution is the coexistence of 
two homogeneous phases whose densities are the smallest 
and largest solutions of Eq. H106|) . 

We now will calculate ec by determining the range of 
e for which the equation G'{p) = has two solutions for 
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p in [0, l/(?]. We have 



G'{p) = 



2 



1 -Pc 



e(2cr- 1) 

+ 1 - Pco: 



dp 



(108) 



By differentiating Eq. H103|) with respect to p one can 
show that dpcon/dp < for pcon in [0, 1] and p in [0, 1/q]. 
Accordingly, setting G"(p) — leads to the quadratic 
equation 

2aep2^„ - pco„(l + 2e + 2ae) + 2e + 2 = . (109) 

This equation has two real roots for either e > 6+ or 
e < e_ and no real roots otherwise, where e± satisfies 



(l + 2e + 2(Te)2 = 16cre(e+l) , 



which gives 



3(7 - 1 ± 2^cr(2cr - 1) 

2(^^^T? 



(110) 



(111) 



Since the two values of pcon corresponding to e_ are both 
greater than 1, we see that two relevant solutions for 
G'(p) = occur only for e > ec = €+. From the value of 
Ec we obtain the transition temperature Tc for the liquid- 
gas transition for interacting dimers from 



1 



(112) 



We close this section by discussing the accuracy of the 
"strong form of tree decoupling" in which we approxi- 
mate a d-dimensional lattice by a Bethe lattice of the 
same coordination number. Recall that in the tree de- 
coupling introduced above Eq. (|11|) it was only assumed 
that we could neglect indirect paths which connect near- 
est neighboring sites A and B. A comparable approxima- 
tion here is to assume that the conditional probability 
that a neighbor of site A is vacant given that both sites 
A and B are vacant is the same as the conditional prob- 
ability when only site A is vacant. If this were the only 
approximation, then the present treatment of the inter- 
acting dimer system would be expected to be quite accu- 
rate - as is the tree decoupling for noninteracting dimers^ 
However, here we also started from the much stronger as- 
sumption that the random variables Xi introduced above 
Eq. H92|l are independent of one another. On the Bethe 
lattice this is obviously true, as can be seen from Fig. ^ 
However, on a real lattice, this assumption neglects the 
many next-nearest neighbor connections between the Xi 
sites. This same neglect reappears if one tries to apply 
the subsequent calculation of p' to a d dimensional lat- 
tice. Accordingly it would not be surprising if the Bethe 
approximation for interacting dimers had an accuracy 
similar to that of mean field theory for the Ising model 
in the same spatial dimensionality. In two dimensions, 
this would imply that the Bethe lattice value of Cc may 
differ from the exact value for a d dimensional lattice by 
about 30%. 



B. Solution by Construction of an Effective 
Hamiltonian 



In this section we generalize the effective Hamiltonian 
for noninteracting dimers to include dimer-dimer interac- 
tions. We continue to mark sites with operators Si which 
obey the trace rules of Eq. I|15|l . But in addition we 
need to have operators which keep track of interactions. 
So at each site i we introduce operators ti which have 
zero trace unless accompanied by an Si operator. These 
operators commute with one another and obey 



(113) 



for p and r each assuming the values 0, 1, ... g. Now we 

set 



n 



(1 -I- xsiSj) (1 + 5titj) 



(114) 



For a configuration of Njj dimer bonds and TV/ interact- 
ing bonds this Hamiltonian gives a contribution to the 
partition function of 



[x(l-f (5)]^°[1-H(5]^^ 



(115) 



To get the desired partition function we thus set 6 = 
e"/3 _ X = g and x = ze"''". 

We now obtain the exact solution for the Bethe lat- 
tice. We could introduce replicas to obtain an expansion 
for the free energy rather than for the partition func- 
tion. However, in the interest of simplicity we work with 
the partition function to obtain an exact solution for the 
Bethe lattice. Accordingly Eq. (|21|l in this case assumes 
the form 



Tr, 



The solution to this equation is of the form 

gi = A + Bs^ + CU + DsiU . 

To satisfy Eq. (|116|) the constants must obey 

A = Tr,5j/Tr,gJ 

B - ze-''"Tr,s,5j/Tr,5« 

C = eTr,i,5j/Tr,5| 



(116) 



(117) 



(118) 



These equations may be written in terms of A, rg 
B/A, rc = C/A and ro = D/A, as 



A' = 



l + q{rB+rD){l + rcY 
l + (7{rB+rD){l+rcY^^ 
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rc 



td = 



l + (j{rB +rD){l+rc) 
From the equation for rc we obtain 



— = erB. (119) 



TB + TD 



(120) 



Using this in conjunction with the equation for tb we 
obtain an equation which determines rc'- 



r-c/(l + e) 



z{l + rcY 



a{e - rc){l + rcY 



(l + e)[l + rc/(e-rc)] 

(121) 



which can be put into the form 

— In (crz) = 2 hi(e — rc) — hi(erc) 

+ (2a- l)hi(l + rc) . (122) 

This is identical to H106|l when we make the identification 
rc/e = Peon- (It has to be admitted that this physical 
interpretation of rc is not obvious if one only has the 
effective Hamiltonian.) 

In Appendix^we develop mean field theory by a suit- 
able decoupling of the effective Hamiltonian and, as ex- 
pected, we obtain results qualitatively similar to those 
for the Bethe lattice. 



IV. QUENCHED RANDOMNESS 

A. Dimers on Percolation Clusters at Infinite 
Fugacity 

Here we consider the statistics of dimers on a quenched 
random lattice in which sites can be either "X" sites with 
probability p or "Y" sites with probability 1 — p. For a 
given configuration C (i. e. for a given distribution of 
X and Y sites), the grand canonical partition function, 
Z{C;{z}) for dimer coverings is calculated as 



Z{C;{z}) = 



Enxx {A) nxY (A) uyy {A) 



'YY 



(123) 



where the sum is over all arrangements ^ of 0, 1, 2, 
. . . hard-core dimers, tiab (-^) is the number of dimers 
covering an A site and a B site in the arrangement A 
(where A and B each assume the values X and Y), and 
ZAB is the activity of an AB dimer. Then the quenched 
average free energy F is calculated as 



F(p;{z}) -^P(C;p) In Z(C;{z}) , 



(124) 



with P(C;p) =p"-(^)(l-p)"«('^), where n^{C) (ny(C)) is 
the number of X (Y) sites in the configuration C. From 



this quenched free energy one then obtains the average 
number of AB dimers as 



Nab(j>-Az\) 



ZAB- 



d F{p-{z\) 

dzAB 



(125) 



An exact result for the dimer density on a Bethe lattice 
does not seem easy to obtain. However we now give an 
exact solution for the dimer density on a Bethe lattice 
in the limit zxx oo and zxy — zyy — > 0. Thus we 
consider Md (p) , the average of the maximal number of 
dimers which can be placed on percolation clusters of X 
sites. Even in this limit the result is not trivial because 
the constraint that dimers do not overlap plays a variable 
role depending on the compactness of the cluster. 

To obtain the exact solution we will explicitly evaluate 
the expansion of M]j(p) in powers of p. For this purpose 
we temporarily consider as a function of the set of 
Pi's, where pi is the probability that site i is present (i. 
e. is an X site). Then we write 



MD{{p^})=Y.^^D{0,■■■0,P^,0,...0) 



Mb(0,...0,k,0,...,0,Pj,0,...0) 



-MoiO, ...0,K, 0, ...0) - MoiO, ...0,Pj, 0, ...0) 

+ ... . (126) 
We write this as 

MDi{p^}) - ^MD(Pz)c + 5I^'^i5(P»'Pj)c 



- ^ AlD{pi,Pj,Pk)c 
i<j<k 



(127) 



where we only indicate as arguments those pi's which are 
nonzero and we introduce the cumulants via 

MD{p^)c = MDip^) 

MD{pi,Pj)c = MD{Pi,P]) - MoiPi) - Md{Pj) , 

(128) 

etc., as in Eq. (47). Note that MoiPi) = be- 
cause a dimer requires two sites being present. Similarly 
MoiPiiPj) 7^ only if sites i and j are nearest neigh- 
bors. On can show that the cumulant vanishes for a 
disconnected diagram. The general term (evaluated for 
Pt = p) is 

= ^A/z5(7)(-l)^--^- 

= The terms of order p^"" in Md{T) 

Nr 



Nr 
— p^ 



(129) 




FIG. 6: Small clusters for which we give the cumulant value 
of M. 

In the first line T denotes a set of pi's, 7 is a subset of 
r (with 7 = r allowed), and Nr is the number of pi's in 
the set r. In the last line Xi — 1 means that the site i 
is included in the set and Xi — —I means that the site 
i is not included. (These definitions follow from the fact 
that to get a term of order p^^, we take a factor of p 
if the site is included and {1 — p) {—p) if the site is 
not included.) For example, for the sets of sites shown in 
Fig. El we find that 

Md{Ti), = p" , 

Md{T2)c = p'[l-2(l)] - V , 

Md{T^)c - /[2- 4(1) + 3(1)]=/ , 

MoiTi), = /[I - 3(1) + 3(1)] = / . (130) 

For r2 there is one term with all the a;'s equal to +1 and 
two nonzero terms with one Xi equal —1. For Fa there is 
one term with all the a;'s equal to +1, four nonzero terms 
with one Xi equal —1, and three nonzero terms with two 
XiS equal to —1. For r4 there is one term with all the x's 
equal to +1, three nonzero terms with one Xi = —1, and 
three nonzero terms with two XiS equal to —1. These 
results suggest that for any connected cluster of sites F 
on the Bethe lattice one has 

Md{T), = {^pf- . (131) 

Note that in contrast to the cumulant, the bare value 
Mo(F) is not simply a function of iVr. 

The proof of Eq. H131|) is by induction on Nr. We 
have explicitly shown this result to be true for A^r equal 
to 2, 3, and 4. We now show that if Eq. (|131|) is assumed 
to hold for iVr < N ~ I, then it holds for Nt = N 
(assuming N > 2). When this is proved, the general 
result is established. 

Consider a diagram with N sites and label the sites so 
that the Nth site is an "end", that is, it is connected to 
only one other site in the diagram and this other site is 
labeled — 1, as shown in Fig. 7. For a Bethe lattice 
(in contrast to the case of hypercubic lattices) this con- 
struction is possible because all diagrams have at least 
one free end. 
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FIG. 7; Cluster of A'' sites in which site A*' is a free end. 
We will use Eq. I|129() for the cumulant, so that 

M_D(rjv)c=P^ ^ XiX2---XnF{xi,X2,-..Xn) , 

(132) 

which we write as 

Md{Tn)c = X! XiX2...Xn-2 

y.\A-B^C -B\ , (133) 

where the algebraic signs reflect the values of the quantity 
{xn-iXn) with 

A = Md{xi,X2 ■ ■ ■XN-2AA) 
B = Md{xi,X2 ■ ■ ■Xiq-2,'^-,-'^) 

C = MD(a;i,a;2 . . .a;Ar_2,-l,-l) 

D = Md{xi,X2...xn-2,-1,1) ■ (134) 

In terms _B, C, and £), the maximal covering by dimers 
does not include the bond — 1 to A^, because for this 
bond to be included obviously both sites A^ — 1 and A^ 
must be occupied (so that xjv~i — xn = 1). Consider 

r(Fjv)=p^ J2 xiX2...XN-2[-B + C] .{135) 

X\ ,X2 . . .2:7V — 2 

This is almost the contribution to the cumulant when the 
site A'^ is not included either in the summations or in the 
covering. We say " almost" because the sum includes the 
factor —pxN which ought to be taken out if we want to 
identify this with the A^ — 1 site problem. We have 

^ X1X2 ■ ■ ■XN-^2[-B + C] 

Xi ,X2...Xn-2 

= {~P)P''-^ Yl XiX2...XN-2[B-C] 

Xi .X2 . . .Xjv — 2 

= {-p)Md{T'j,_^), , (136) 

where r^_-^ is the diagram obtained from Fjv by omitting 
the iVth site. 

Now consider A. Suppose the maximal covering does 
not actually include the bond A^ — 1 to A^. Then clearly, 
if this covering is to be maximal it must include a bond 
from some site to A^ — 1. But now we may keep the 
number of dimers maximal by moving this bond which 
includes the site A^ — 1 to cover the bond A^ — 1 to A^. 
So, without loss of generality, in the term A the maximal 
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covering can be chosen so as to include the bond N — 1 
to N. Thus for > 1, 



„N 



P 



X1X2 ■ ■ ■ XM-2A 



X1X2 ■ ■ ■ Xn~2 

x[Md{xi,X2,...xn-2) + 1] ■ (137) 

Now if iV > 2 the term with the 1 wiU vanish when 
summed over xi. (For N — 2 there are no summations 
left and this is a special case.) Then we see that for 
N >2 



P 



E 

Xl ,X2---XN — 2 



XIX2..-XN-2A = p'^Md{T%_2)c , 



(138) 

where F'^ ^ is obtained from F^v by deleting sites N and 
A'^ — 1. The same reasoning can be applied to term D, 
but in this case there is no term with 1. So 



P 



^XiX2 . ■ ■XN-2{-D) = -p^Af_D(F^_2)c , 



(139) 



Combining all our results we have (for > 2) 



MDiTN)c = P^ J2 XiX2...Xn[A-B + C-D] 



Xl ,X2 ...Xjsf — 2 



(-p)Afz5(F^_i)c 



(140) 



This completes the proof by induction. 

For a general lattice we may make the tree approxima- 
tion in which the sum is carried over all diagrams with no 
loops. In that case the result is that the average dimer 
density is given by 

M{p) = AId{p)/No = N^'Y.(-p)''" ' (141) 



where A'o is the total number of sites, the sum is over 
all loopless clusters C of 2 or more sites, and Nc is the 
number of sites in that cluster. Thus for the Bethe lattice 
Eq. ((TTTjl yields 



(142) 



n>l 



where W{n) is the number (per site) of clusters of n sites 
on a Bethe lattice. (To avoid edge effects it may be more 
precise to say that nW{n) is the number of clusters of n 
sites, one site of which is the central site in an arbitrarily 
large Cayley tree. This definition ensures that this result 
applies to hypercubic lattices in the asymptotic limit of 
large dimensionality.) We use the result of Fisher and 
Essami^ 



W{n+l) 



(<j + l){an + a)\ 



(n + l)!((Tn - 



1)! 



(143) 



Therefore the exact result can be written as 

(a + l)(an)! ^ 
n\(an — n + 2)! 



(144) 



n>l 



This result indicates that the singularity in M{p) is at 
p — —pc, where Pc — {(J — /(j"^ is the critical con- 
centration for branched polymers on the Bethe latticei2£ 
Curiously then, dimer statistics on percolation clusters 
is related to a somewhat artificial model of localization 
(which is also related to branched polymer statistics in 
the same wayi^) 

For (T = 1 (a linear chain) this gives 



Mip) 



P 



1+p 



(145) 



This result could be obtained far more simply by noting 
that the expected number of clusters per site of length n 
sites is p"(l — p)^ and for a cluster of length 2n or 2n+ 1 
the maximal number of dimers is n, so that 

00 

M{p) = (P'" i^-pfn , (146) 

n=l 

which reproduces the result of Eq. H145|) . 



B. Alternate Expression 

We now give a closed form expression for M{p). To 
do that we construct an expression for cPM/dp^. In this 
quantity the binomial coefficient can be expressed as a 
contour integral: 



(a + l)£(-p)' 



dz 1 + z)'"" , 



where the contour surrounds the origin. The sum over n 
is a geometric series which can be summed. After some 
manipulations we found that 



M{p) = p + zo 



-Zq , 



(148) 



where zq is the root of the algebraic equation 

(1 + zoY = -zo/p , (149) 

which is proportional to p for small p. One can check 
that this gives the correct result for = 1. For ct = 2 it 
gives 



Mip) 



1 

4p2 



1 - 6p - 6p^ 4- 4p-^ + + 4p) 



^3/2 



(150) 



This result shows that there is no singularity on the pos- 
itive p axis. 
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RANDOM DEPOSITION 



Introduction 



Here we consider a special case of a model of catalysis^ 
in which hard-core dimers are randomly deposited on 
bonds and the hard-core constraint does not allow two 
dimers to intersect the same site. Deposition on all al- 
lowable bonds is equiprobable. A quantity of interest is 
the final concentration when no further dimers can be 
deposited. We have not been able to construct such a 
solution on a Bethe lattice. However, here we give an 
exact solution of this model in one dimension. 

To characterize this process it is essential to consider 
deposition on a line of bonds of finite length in one di- 
mension. In order to consider deposition recursively, we 
therefore introduce the function F{i,j), (with i < j) 
which is defined to be the average number of dimers 
which will ultimately be deposited in the interval between 
two preexisting dimers, one on bond i and the other on 
bond j. Obviously F{i, j) = F{i — n,j — n) and as short- 
hand we set F{l,j) = F{j). 

Let us see what this function is for small argument. It 
is clear that F(l,3) = F(l,4) = because a new dimer 
can not be deposited on a bond neighboring an occupied 
bond, since neighboring bonds share a site which can 
not be occupied by two hard-core dimers. Next consider 
F{1, 5) = F{5). If we start with dimers on bonds #1 and 
#5, then an additional dimer can only be deposited on 
bond #3, so 



F(1,5) = F(5) = 1 



(151) 



is the average number of dimers which will be deposited 
between dimers at sites #1 and #5. Similar considera- 
tions indicate that F{6) — 1 and F{7) — 5/3. 

Next consider F(l,8). In the first step an additional 
dimer will be placed on bonds #3, #4, #5, or #6, each 
with probability 1/4. If it is placed on bond #3 (or 
equivalently on bond #6) , then we have added one dimer 
and will be able to add (on average) F(l, 3) -I- F(3, 8) = 
+ 1 = 1 further dimer in later step(s). Therefore in each 
case these two processes lead to the deposition of two 
dimers. So the combined contribution to F(l,8) from 
these two cases, each occurring with probability (1/4) is 
dF(l,8) — 1. Similarly, if the first new dimer is placed 
on bond #4 (or equivalently on bond #5), then we have 
added one dimer and will be able to add (on average) 
F(l,4) -I- -F(4, 8) = -f 1 = 1 further dimer. As before, 
the contribution to F{1, 8) from these two cases, each oc- 
curring with probability (1/4) is 6F{1,8) = 1. Therefore 
we see that F(8) = 2. 



B. Recursion Relation 

Now imagine starting with dimers on bonds #1 and , 
with > 4. The first added dimer can go on bond #3, 



#4, .... - 3, - 2, each with probability 1/(A^ - 4). 
So, if we include this added dimer we have 

Fil.N) = l + -r— - J2 m,j) + F{j,N)] 



= 1 



N -4 
2 



J=3 



N -A 



(152) 



J=3 



Also 



F{N - 1) = 1 



Using these we see that 



A^- 5 



E ^C?)- (153) 



J=3 



{N-4)F{N) = (A^ - 5)F(A^ - 1) + 1 

+2F{N-2). (154) 



Now we form the generating function 

oo 

F(x) = ^F(A^) 



N 



(155) 



^=5 



Now multiply Eq. (|154|l by and sum from A^ = 5 to 
N ^ oo. Keeping in mind that F{N) vanishes for A^ < 5 
we get 



xF^ -AF = x^F^ - AxF - 



1-x 

where F^ = dF(x)/dx. From this we find that 



2x^F , (156) 



Fix) 



(157) 



{l~xf\ 2 J ' 
To see what this implies about F{N) for large A^ we write 

^ ' (1 - .X)2 

^ Hjl) gi(l) 
(l-a;)2 1-x 

(158) 

where i?„(l) = (PH{x)/dx"\x=i- From this we see that 
F(A^) - i/(l)Ar = (1 -(5)(A^/2) , (159) 

where 

5 = « 0.1354 (160) 

is the fraction of sites which remain vacant in the jam- 
ming limit (after deposition is completed). 
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VI. SUMMARY 

In this paper we have presented two approaches, which 
we call geometric and algebraic, for the analysis of dimer 
statistics. Both approaches yield exact results when the 
lattice does not have loops, as for a Bethe lattice. The 
charm of the geometrical method is that it starts from the 
most basic statement, namely that in a grand canonical 
ensemble, the activity z is precisely the ratio of the statis- 
tical weights of the TV-particle system to that of the A^-l-l- 
particle system. The virtue of the algebraic approach is 
that once the Hamiltonian is constructed using operators 
whose trace rules incorporate kinematic restrictions, the 
standard procedures of statistical mechanics can be ap- 
plied, for instance to obtain series expansions for finite 
dimensional lattices. (This approach was used^^ to give a 
field theoretic analysis of the monomer-dimer problem.) 
Here we also show that the use of replicas (which nor- 
mally are invoked to implement quenched averages of a 
random Hamiltonian) can be useful in converting a series 
for the partition function into one for its logarithm, the 
free energy. Exact solutions for generalizations in which 
a) the lattice is anisotropic or b) dimer-dimer interactions 
are included were also developed using both approaches. 
We give an exact solution for dimer statistics on a Bethe 
lattice in a simple quenched random potential in which 
dimers are placed on percolation clusters. Finally, we de- 
veloped an exact solution for the fraction of sites which 
remain vacant after random disposition is completed in 
a one dimensional system. 
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sponding density of dimers pi, p2, and pa, then we have 
ln[P(p3)/P(pi)] = ^ \G{p) - ln(l/a^)] dp .(A2) 

If zo (e) is the value of z such that the area of the loop 
below the dashed line at height ln(l/(Tz) is equal to the 
area of the loop above that line, then P{p^)/P{pi) = 1. 
If z is slightly less than zq , then the negative area exceeds 
the positive area and P{p^) / P{pi) < 1, in which case is 
the stable phase has p = pi. Note that p2 is always 
less probable than pi. Similarly, if z is slightly greater 
than Zq z = z<, then the stable (most probable) phase 
corresponds to p = p^. Note that because of the pref- 
actor Nq/2 before the integral, "most probably" means 
"overwhelmingly most probable." 

If we place a density of dimers on the lattice which is 
intermediate between pi{zo) and Ps^zq), the lattice will 
separate into regions with coverage at densities pi(zo) 
and ps^zq), the size of the two regions being determined 
by the requirement that the total number of dimers is 
that specified. 

If e is slightly greater than Ec, then the temperature T 
is slightly less than Tc and the definition of e implies that 
(Tc — r) cx (e — Ec). In this case the distance between the 
two roots of Eq. (|109|l is proportional to (e — Ec)^/^, and 
thus the distance between the corresponding value of p 
is also proportional to (e — ec)^/^, which is proportional 
to Tc — r)^/^. Examination of Fig. |S1 without additional 
analysis, makes it clear that the difference in density be- 
tween the two coexisting phases, 2/3(29) — yi(zo) is also 
proportional to {Tc - Tf/'^. 

APPENDIX B: MEAN FIELD THEORY 



APPENDIX A: ANALYSIS OF THE PHASE 
TRANSITION 



The phase transition in the presence of dimer-dimer in- 
teractions is entirely similar to that in the van der Waals 
gas. 

The right-hand side of Eq. (|105() gives an explicit for- 
mula for zP{Nd + 1)/P{Nd) when Eq. fTU^ is usedM. 
Thus we have 

HPiNo + l)/PiND)] = G{p) + In(az) . (Al) 

Figure [SI shows schematic graphs of G{p) versus p for e 
slightly less than and slightly greater than tc- For a = 2 
we find Cc « 4.95 and the value of p at the critical point 
is pc « 0.159, G{pc) = 3.226, and Zc « 0.0199. 

If e > ec and z is close to Zc, a horizontal line (dashed 
in Fig. O at height ln(l/(Tz) will intersect the graph of 
G{p) three times, each of the intersections corresponding 
to a value of No for which P{Nd + 1)IP{Nd) = 1. If 
we denote the intersections yi < 1/2 < 2/3 and the corre- 



In this appendix we obtain mean field theory for inter- 
acting dimers from the Hamiltonian using the standard 
decoupling even though this Hamiltonian involves oper- 
ators which obey unusual trace rules. 

In view of the trace rules (extended to infinite q) we 
may write 



,zsiSj-{-af3titj 



(Bl) 



where z = ze Mean field theory is obtained by ig- 
noring correlated fiuctuations and writing 



(3n = J2 



z(s,(s) +Sj{s) - (sf) 



Then 



(5) 



+af3 {t,{t)+tj{t) - (t)2) 



^-^^g^QSiqz{s)+qal3{t)ti 
^^^8,qz{s)+qat}{t)ti 



(B2) 
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1 + qz{s)ei"l^<*^ 



(B3) 



D 



(t) 



^-^^j-^^Siqz{s)+qaf3{t)ti 
'J'y..QSiqz{s)+qaf3{t)ti 

1 + qz{s)ei°'P(*') 



B 



We solve Eq. IIB4p for (s) to get 



(B4) 



(B5) 



so that Eq. (|B3p gives 

\n{t) ^ \n{qz) + 2qa/3(t) + 21n[l - {t)] . (B6) 

To shorten the discussion we assume a second order 
transition in which case the above equations and its two 
derivatives with respect to {t) are zero. Differentiating 
twice we get 



1 



l-(t) 



and 



(B7) 



(B8) 



The last equation gives tc, the value of {t) at the critical 
point, to be tc = \/2— 1- Putting this into the preceding 
equation gives qajSc = + (3/2). 

We may compare this with the solution for the Bethe 
lattice. For simple models (such as the Ising model) it 
is known that the Bethe lattice solution only coincides 
with mean field theory for large q. For large q the Bethe 
lattice solution for Tc given by Eq. (|112|l agrees with the 
present result. 



APPENDIX C: GEOMETRICAL CALCULATION 
OF FIRST LOOP CORRECTION FOR THE 
TRIANGULAR LATTICE 

For the triangular lattice, a simple geometrical argu- 
ment permits us to calculate most of the difference be- 
tween the "exact" value of N o/N (notation as in Table 
I) and the value given by the Bethe approximation. The 
argument depends on the "tight" structure of the trian- 
gular lattice, and we have not been able to extend the 
argument to the square lattice. Furthermore, unlike the 
"loop corrections" calculated in Sec. IIC, our geometrical 
calculation is not the first term in a systematic series of 
corrections. Nevertheless, the calculation is simple and 
remarkably accurate. From Eq. (7) we have 



1/^ = Pv/p , 



(CI) 



G 



FIG. 8: A portion of the triangular lattice. 

where pv is the probability that adjacent sites A and B 
are both vacant and p is the probability that a bond is 
occupied. Clearly 

P(B vacant) = P(B vacant and A vacant) 

+P(B vacant and A occupied) (.C2) 



Consider the section of the triangular lattice shown in 
Fig. IHI We see that 

P(B vacant) = P(A vacant) = 1 - 6p (C3) 

and we set 

P(B vacant and A vacant) = (1 - 6p)(p') , (C4) 
where 

p' = P(B vacant | A vacant) . (C5) 

Furthermore 

P(B vacant and A occupied) 

= P(B vac and AE occ) -f P(B vac and AD occ) 
+P(B vac and AF occ) + P(B vac and AC occ) 
-f P(B vac and AG occ) . (C6) 

We now estimate each of the terms on the right-hand 
side of this equation. It seems evident that there 
is very little difference between P(B vacant | AE occ) 
and P(B vac| A vacant), since the influence of the ex- 
tra dimer would have to propagate over a long tortuous 
path. Thus we write P(B vac and AE occ) = {p){p'). 
Similarly, and on slightly weaker ground, we write 

P(B vacant and AD occupied) 
= P(B vacant and AF occupied) = (p)(p') . (C7) 

To estimate P(B vac and AC occ) we note that if AC is 
occupied then AB and CB must be unoccupied. As far 
as state of site B is concerned, it makes little difference 
whether we specify that bond AC is occupied or that 
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bonds AB and CB are unoccupied, because the effect of 
the presence of the dimer AC on the state of site B is 
very indirect. Thus we make the approximation that 

P(B vacant | AC occupied) 
= P{B vacant | AB and BC unoccupied) . (C8) 

However, the conditional probabihty that B is vacant, 
given that AB and CB are unoccupied, is (1— 6p)/(l— 2/7), 
so that approximately 

P(B vac I AC occ) = (1 - 6/9)/(l - 2p) , (C9) 

and thereby that 

P(B vac and AC occ) = p(l - 6p)/(l - 2p) (CIO) 

Now we use our evaluation of each term on the right-hand 
side of Eq. ljC6|l to write Eq. ljC2|l as 

(l-6p) = {1 - 6p){p') + {3p){p') 



+ (2p)(l-6p)/l-2p) . (Cll) 

Solving for p' we find that 

p' = (l-6p)(l-4p)/(l-2p)(l-3p) (C12) 

and Eq. yields 

1/z = (l-6p)(p')/(p) 

= (1 - 6p)2(l - 4p)/(p){l - 2p){l - 3p) (C13) 

The computer readily calculates p for each value of z 
in Table II. The corresponding value of N u/N , which we 
call ttq, is 3p and is exhibited in Table II. 

Thus, a simple geometrical argument yields a good es- 
timate of the "first loop correction" for the triangular 
lattice. However, we have been unable to make a cor- 
responding argument for the square lattice, nor can we 
extend this argument in an orderly way to make higher 
order corrections. 
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